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Abstract 

A model for Galactic chemical evolution, driven by supernova-induced star formation, is formulated 
and used to examine the nature of the Galactic halo at early epochs. In this model, new stars are formed 
following each supernova event, thus their abundance pattern is determined by the combination of heavy 
elements ejected from the supernova itself and those elements which are already present in the interstellar 
gas swept up by the supernova remnant. The end result is a prediction of large scatter in the abundance 
ratios among low-metallicity stars, reflecting a different nucleosynthesis yield for each Type II supernova 
with a different progenitor mass. Formation of new stars is terminated when supernova remnants sweep up 
too little gas to form shells. We show from calculations based on the above scenario that (i) the observed 
[Fc/H] distribution for the Galactic halo field stars can be reproduced without effectively decreasing the 
heavy-element yields from Type II supernovae by some manipulation required by previous models (e.g., via 
mass loss from the early Galaxy, or later mixing with "pristine" hydrogen clouds), (ii) the large observed 
scatter in the abundance ratio [Eu/Fe] for the most metal-poor stars can also be reproduced, and (iii) the 
frequency distribution of stars in the [Eu/Fe] — [Fe/H] plane can be predicted. Our model suggests that 
the probability of identifying essentially metal-free stars (Population III) in the local halo is around one in 
10 3-4 , provided that star formation in the halo is confined to individual gas clouds with mass of 10 6-7 Mq 
and that the initial mass function of metal-free stars is not significantly different from the Salpeter mass 
function. 
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1. INTRODUCTION 

In conventional chemical evolution models, stars are as- 
sumed to form from well-mixed gas clouds at a rate pro- 
portional to some power of the gas density (Schmidt 
1959), thus inheriting the abundance pattern of the gas 
at that time (e.g., Tinslcy 1980). This simplified treat- 
ment of the complex star formation process has been re- 
markably successful in accounting for the general features 
of the chemical compositions of nearby stars and H n 
regions in a consistent manner (e.g., Pagel & Patchett 
1975; Matteucci & Greggio 1986; Yoshii, Tsujimoto, & 



Nomoto 1996). However, Shigeyama & Tsujimoto (1998, 
hereafter ST98) and Tsujimoto & Shigeyama (1998, 
hereafter TS98) argued that these conventional models 
cannot be applied to the early stage of the Galactic 
halo, because observed abundance patterns of extremely 
metal-deficient stars in the range of —4 <[Fe/H]< —2.5 
(McWilliam et al. 1995; Ryan, Norris, & Beers 1996) are 
incompatible with those predicted by these simple for- 
mulations, (see also Audouze & Silk 1995). 

As an alternative, they proposed a scenario in which 
stars are born from each supernova remnant (SNR), fol- 
lowing a sweeping up of the interstellar gas and formation 
of a dense shell immediately behind the radiative shock 
front. The dense shell is dynamically overstable (Vish- 
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niac 1983; Ryu & Vishniac 1987), so that it is broken 
into a few thousand fragments, in each of which the self- 
gravity is unimportant. According to Nakano (1998) it 
is likely that, for at least some of these fragments, the 
ambient pressure outside the fragments exceeds the crit- 
ical value above which no equilibrium state exists. Such 
fragments will dynamically contract and eventually form 
stars. 

This hypothesis sheds light on the entire process of 
early-epoch star formation. Galactic chemical evolution 
is envisioned as a successive sequence of supernova (SN) 
explosions, shell formation, and resulting star formation. 
In this picture, the abundance pattern of each star is set 
not only by the heavy elements ejected from the SN ex- 
plosion which directly preceeds local star formation, but 
also incorporates those elements already present in the 
interstellar gas which are swept up by the SNR. The pre- 
dicted stellar abundance patterns are thus different from 
those of the gas at the time when stars form. This dif- 
ference can be quite large in the early stage of Galactic 
evolution, when the metallicity in the gas is very low. 
The abundance ratios of low-metallicity stars are pre- 
dicted to exhibit a large star-to-star scatter, depending 
in detail on the abundance patterns of SN ejecta with 
different progenitor masses. 

By incorporation of a sequence of SN-induced star for- 
mation into the enrichment process of heavy elements, we 
describe the evolution of stellar abundances in the Galac- 
tic halo characterized by the inhomogeneous mixing of 
chemical compositions in dense shells and the ambient 
medium. 



2. CHEMICAL EVOLUTION IN A CLOUD 

We assume that the star-forming process is confined to 
separate clouds, of mass M c , which make up the entire 
halo. In this section we present a formulation to describe 
the chemical evolution in a cloud that is initially com- 
posed of metal-free Population III stars (Pop III) and 
gas which has yet to form stars. The mass fraction of 
such stars in our model is a parameter hereafter denoted 
by xiu, the massive ones of which eventually explode as 
Pop III SNe to initiate chemical evolution. 

All stars of subsequent generations are assumed to 
form in SNR shells. The mass fraction of each shell which 
turns into stars is assumed to be constant, and is denoted 
by e. Heavy elements ejected from an SN are assumed to 
be trapped and well-mixed within the SNR shell. Some 
of these elements go into stars of the next generation, 
and the rest is left in the gas that will be mixed with the 
ambient medium. The above process will repeat with in- 
creasing metallicity until SNRs can no longer sweep up 
enough gas to form shells. No stars will form from SNRs 
after this happens, and the process terminates. 



When stars form as above, the star formation rate 
(SFR) M*(t) induced by Pop III SNe at time t measured 
from the birth of their progenitor stars is given by 

drat 



M*(t) = eM sh (m t ,t) 



dt 



x m M c , 



(1) 



and the SFR for later generations can be expressed as 

f mu Mm) ■ 

M*(t)= dm€M sYl {m,t)-^-M- t {t-T{m)), (2) 

imsx(ro t ,msN,i) m 

where r(m) denotes the lifetime of a star of mass m, 
and m t is the stellar mass for which r(m) = t. The 
initial mass function (IMF) <p(m) used here is a Salpeter 
one, with a slope index of —1.35. An upper mass limit 
of stars is assumed to be m u — 50 M Q ; a lower mass 
limit of mi — 0.05 M Q (Tsujimoto et al. 1997) is also 
assumed. The lower mass limit for stars that explode as 
SNe is taken to be tosn.z = 10 M Q . The mass of the shell 
M s h (to, t) formed at time t from an SN with progenitor 
mass m is a sum of the mass of the SN ejecta M e j(m) 
and the mass of the swept-up gas M sw (m, t); 



M sh (m,t) = M ej (m) + M sw (m,t), 
with 

M sw (to, t) = p g (t)R 3 SN (E (to) , p g (t)) , 



(3) 



(4) 



where p g {t) is the density of the interstellar gas and 
-Rsn (E(m), p g {t)) is the maximum radius of the SNR 
shell. M sw (m,t) is insensitive to time t because M sw 
is proportional to p g (t)-° m (see ST98). The SN ex- 
plosion energy is assumed to be E{m) = 10 51 ergs, ir- 
respective of stellar mass to. Thus a constant value of 
M sw (m,t) — 6.5 x 10 4 M Q will be used below. 

Utilizing the SFR introduced in equations (1) and (2), 
the mass M g {t) and the metallicity Z g (t) of the gas will 
change with time according to the formulae: 



dM g 
dt 

and 



Mm) 



= -M*{t) + / dmM cj (m)^^M*(t - t(to)),(5) 
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d(z M ) f mu 

— JT^ = -j ^.(m,t)eM s h(m,i) 

Jmax(mt.msNj) 



(j)(m) 



TO 



<p(m) 
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x dm'eM s h(m' ,t — r(m)) 



0(m') 



to' 



M*(i-r(m) -r(m')),(6) 



where M z (m) is the mass of all synthesized heavy ele- 
ments ejected from a star with mass m, and Z*(m 7 t) is 
the metallicity of stars born at time t from an SNR shell 
with progenitor mass to, which is defined as 

M z (w) +Z g (t)M sw (m,t) 
M cj (m) + M sw (m,t) 

Since M*(0) = x m M c and M g (0) = (1 - x m )M c at 
t = 0, after specifying xm and M c and setting Z g (0) = 
0, equations (5) and (6) are then integrated with (1)- 
(4), and (7). Star formation stops when the total mass 
involved in SNR shells reaches the total mass of gas, i.e., 



Z*(m, t) 



(7) 



J t—AtJ max(m ( / , mg N . ;) 



m 



r(m)) = M s (t), (8) 



where At is the time for the turbulent motion in the 
ambient medium to completely destroy the shape of an 
SNR shell. When this condition is satisfied, the cloud is 
composed of stars, hot tenuous gas with negligible mass, 
and dense blobs of gas (an SNR shell would be broken 
into blobs) which fill in a negligible volume of the cloud. 
In the following calculation At = 3 x 10 6 yr is adopted 
which corresponds to the crossing time of the turbulent 
flow over an SNR. 

The heavy-element yields from SNe are estimated di- 
rectly from the observed abundances in metal-poor stars, 
as opposed to the presently uncertain predictions of theo- 
retical supernovae models, using the procedure proposed 
by ST98 and TS98. We note that the europium yield is 
estimated from the data of Ryan et al. (1996). 



3. RESULTS 

From the formulation in the preceding section, we cal- 
culate the chemical evolution for Galactic halo stars. The 
free parameters in our model are the mass fraction im 
of metal-free Pop III stars initially formed in each cloud 
and the mass fraction e of stars formed in a dense shell 
of each SNR. The values of these parameters are chosen 
to reproduce the observed [Fe/H] distribution function 
of halo field stars for [Fe/H]< —1. If ini is larger than 
about 1 x 1CP 2 , the total gas swept up by the first SNRs 
exceeds the entire amount of available gas, and thus star 
formation stops at the first or second generation. In order 
for the enrichment process to continue as a sequence of 
SN-induced star formation, xm must be less than about 
1 x 1CP 2 . In addition, e must be strictly confined within 
a narrow range so that nearly one (slightly larger than 
one) massive star is born from each SNR. This might im- 
ply that star formation is suppressed and self-regulated 
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Fig. 1. The relation between the duration AT of star 
formation in a cloud and a probability pni of ob- 
serving Pop III stars as a fraction of their num- 
ber among all the long-lived stars that have ever 
formed with m < IMq. The dashed lines denote 
the strict lower bound of pni for each mass M c of 
a cloud. 



by the photoionization from an OB star. A value of 
e = 4.3 x 10~ 3 is found to give the best fit to the ob- 
served [Fe/H] distribution function for various values of 
xui(< 1 x 1CP 2 ). If e slightly exceeds this critical value, 
star formation will soon stop with little enrichment by 
heavy elements. On the other hand, if e is taken slightly 
below its critical value, star formation will proceed until 
most of the gas is used up. In such a scenario, the major- 
ity of halo stars become too metal-rich to be consistent 
with observations. 

The duration AT of star formation in a cloud is related 
to xm in the sense that a smaller xm gives a longer AT. 
Here we define a probability pm of observing Pop III stars 
as a fraction of their number among all of the long-lived 
stars that have ever formed with to < IMq. The prob- 
ability pm decreases with increasing AT - this relation 
is shown by the thick line in Figure 1. Each value of pm 
corresponds to about ten times xm, because about 10 % 
of -Mg(O) has been converted to stars in the end. We note 
that the most metal-poor halo stars exhibit abundance 
patterns of genuine SN II origin, showing the overabun- 
dance of a-elements relative to iron Fe (e.g., Wheeler, 
Sneden, & Truran 1989), which means that they were 
born before SNe la start to contaminate the gas with 
their products. If the lower bound for the lifetime of SN 
la is set at 1 Gyr (Yoshii, Tsujimoto, & Kawara 1998), 
the duration AT has an upper bound, and the probabil- 
ity pm should be larger than 10~ 4 . However, such an 
argument does not hold if relatively few SNe la occur at 
[Fe/H]< -1, as argued by Kobayashi et al. (1998). 

One absolute condition under which the SN-induced 
star formation proceeds in a cloud is that at least one 
massive star of Pop III must initially explode as an SN. 
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Fig. 2. The frequency distribution of the long-lived 
halo field stars for [Fe/H] < — 1 against the iron 
abundance, compared with the observations of 
Ryan & Norris (1991). The observational data for 
[Fc/H]> —1 are not shown because these data do 
not correspond to the number pure halo stars due 
to an enormous contamination by disk stars (see 
Ryan & Norris 1991). The model curves have been 
convolved with Gaussians having <r=0.3 dex (solid 
line) and (7=0.15 dex (dotted line), respectively. 



Therefore, for each mass M c of a cloud a lower bound 
of pni exists, and is indicated in Figure 1 by the verti- 
cal dashed line. If M c is in a range of several 10 6_7 M Q , 
which is inferred from the mass scale of the Galactic glob- 
ular clusters, one Pop III star might be observed in a 
complete sample of 10 3 ~ 4 halo stars. This estimate de- 
pends critically on the assumed IMF for Pop III stars. 
For instance, if we adopt the theoretical IMF for metal- 
free stars proposed by Yoshii & Saio (1986), the expected 
piu is reduced by a factor of 15 — 45. 

Figure 2 shows the predicted stellar [Fe/H] distribu- 
tion function for un = 1CP 6 , as compared with the data 
obtained by Ryan & Norris (1991). The result has been 
convolved with a Gaussian of a = 0.15 dex (dotted line) 
and a = 0.3 dex (solid line) because measurement er- 
rors in [Fe/H] lie between these values (Ryan & Norris 
1991). Rather good agreement with the data is obtained 
for [Fe/H]< —1 using the Salpeter IMF, without decreas- 
ing the heavy-element yield by some manipulation such 
as mass loss from the halo (Hartwick 1976; Searle & Zinn 
1978; Bond 1981; Laird ct al. 1988; Ryan & Norris 1991). 
The termination of star formation at the gas metallicity 
of [Fe/H]^ —1.5 produces the peak of stellar frequency 
at [Fe/H]~ —1.6 as observed. Only about 10% of the 
initial gas has been converted to halo stars, and the rest 
is left as the remaining gas which may fall onto the disk 
to be consumed in formation of disk stars. 

Figure 3 shows the predicted [Fc/H] distribution func- 
tion of stars formed at a given age in the age— [Fe/H] 
plane. Three distribution functions shown at different 
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Fig. 3. The predicted [Fc/H] distribution functions of 
the long-lived stars at different ages of 4 X 10 8 , 8 X 
10 s , 1.2 X 10 9 yrs, in the age-[Fe/H] plane. The 
total number of stars in each case is normalized to 
100. The dashed line denotes the age-metallicity 
relation of the gas. 



ages of 4 x 10 8 , 8 x 10 8 , 1.2 x 10 9 yrs, normalized in 
such a way that the total number of stars is 100. The 
dashed line denotes the age-metallicity relation of the 
gas. It is important to note from this figure that the 
stellar metallicity does not correspond to a unique age, 
but that the scatter in stellar [Fc/H] among stars formed 
at a given age progressively diminishes with time. As can 
be appreciated from inspection of the figure, some stars 
with abundance [Fe/H]> —3.0 will form at essentially the 
same time as stars with abundances much lower than this 
value. 

Our primary concern is whether the observed stellar el- 
emental abundance patterns in the most metal-deficient 
halo stars can be predicted by our model. The abun- 
dance ratio [Eu/Fe] is a suitable probe for this pur- 
pose, because this ratio in extremely metal-deficient stars 
(McWilliam et al. 1995; Ryan et al. 1996) spans over 
2 dex, far exceeding the measurement errors. Figure 4 
shows the color-coded frequency distribution of stars in 
the [Eu/Fe] — [Fe/H] plane, normalized to unity when in- 
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Fig. 4. The color-coded frequency distribution of 
the long-lived stars in the [Eu/Fe] — [Fe/H] plane 
convolved with a Gaussian having (X=0.2 dex for 
[Eu/Fe] and ct=0.15 dex for [Fe/H], while there is 
no convolution for the inset. The dashed line in 
the inset denotes the [Eu/Fe] — [Fe/H] relation for 
the gas. The symbols represent the data taken 
from various references (open squares: McWilliam 
ct al. 1995; filled squares: Ryan et al. 1996; plus 
signs: Luck & Bond 1985; open triangles: Gilroy 
et al. 1988; crosses: Magain 1989). 

tegrated over the entire area (see the color bar for the 
scale) . In order to compare with the data, the frequency- 
distribution has been convolved with a Gaussian with 
a = 0.2 dex for [Eu/Fe] and a = 0.15 dex for [Fe/H]. 
For illustrative purposes the frequency distribution with- 
out convolution is also shown in the inset. The pre- 
dicted [Eu/Fe] — [Fe/H] relation for the gas is shown by 
the dashed line. The predicted scatter in [Eu/Fe] be- 
comes smaller toward larger [Fe/H] and converges to a 
plateau at —2 <[Fe/H]< —1. This tendency, which is also 
seen in other elements, is understood by considering that 
the major contribution to stellar metallicity has switched 
from the ejecta of SNe to the metallicity in the interstel- 
lar gas swept up by SNRs (see eq.[7]) in this abundance 
interval. 



a series of star formation events. (3) Star formation ter- 
minates when dense shells cannot be formed from SNRs. 
(4) About 90 % of the cloud mass remains unused in star 
formation and may fall onto the still-forming Galactic 
disk. Our model predicts that stellar abundance pat- 
terns are different from those of the interstellar gas at 
the time when stars form. Calculations presented here 
have shown that the observed features of a large scatter 
in the europium to iron abundance ratios and the metal- 
licity distribution function with ([Fe/H]} ~ —1.6 for halo 
field stars are naturally reproduced. A particularly at- 
tractive feature of our model is that agreement with this 
mean stellar metallicity is achieved without decreasing 
the value of heavy-element yields from SNe. 

Our model strongly suggests the existence of metal- 
free Pop III stars, because it is the first SN explosions 
among them that trigger a series of SN-induced star for- 
mation. Provided that star formation is confined to sep- 
arate clouds with mass of 10 6_7 M Q , and that metal-free 
stars have a Salpeter IMF, it is estimated from the model 
that one Pop III star would be found in a sample of 10 3-4 
halo stars. Ongoing surveys which seek to increase the 
numbers of stars known at the lowest metal abundances 
(e.g., Beers, Preston, & Shectman 1992; Christlieb et 
al. 1999) provide the best hope for the eventual detection 
of this very first generation of stars. However, the proba- 
bility of detecting them would be smaller if the primordial 
IMF is enhanced on the massive part compared to the 
Salpeter IMF (Yoshii & Saio 1986) and/or if the accretion 
of interstellar gas contaminates the surface of metal-free 
stars (Yoshii 1981). Both of these possible complications 
are suitable for direct tests via high-resolution observa- 
tions of a sufficiently large sample of dwarf and giant stars 
with [Fe/H]< —2.0, which already exists (see Beers 1999 
for a review). For an insight into the star formation pro- 
cess, it may be interesting to consider that stars can be 
born from interstellar gas disturbed by shock waves, im- 
plying that the first generation stars (Pop III stars) could 
have been born in a dense shell formed by shock waves 
triggered by a cloud-cloud collision in a protogalaxy. 

This work has been partially supported by COE re- 
search (07CE2002) of the Ministry of Education, Science, 
Culture, and Sports in Japan. We would like to thank 
T. C. Beers for many fruitful discussions. 



4. CONCLUSION 

We have presented a new model for Galactic chemical 
evolution based on the SN-induced star formation hy- 
pothesis. An overview of our model applied to the Galac- 
tic halo is as follows - (1) The metal- free Pop III stars 
form by some (as yet unspecified) mechanism in primor- 
dial gas clouds of the Galactic halo. (2) The most massive 
stars among them explode as Type II SNe, which trigger 
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